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NGC 2579 and the carbon and oxygen abundance 
gradients beyond the solar circle"^ 
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ABSTRACT 

We present deep echelle spectrophotometry of the Galactic Hii region NGC 2579. 
The data have been taken with the Very Large Telescope Ultraviolet- Visual Echelle 
Spectrograph in the 3550-10400 A range. This object, which has been largely ne- 
glected, shows however a rather high surface brightness, a high ionization degree and 
is located at a galactocentric distance of 12.4 ± 0.7 kpc. Therefore, NGC 2579 is an 
excellent probe for studying the behaviour of the gas phase radial abundance gradients 
in the outer disc of the Milky Way. We derive the physical conditions of the nebula 
using several emission line-intensity ratios as well as the abundances of several ionic 
species from the intensity of coUisionally excited lines. We also determine the ionic 
abundances of C^"*", 0+ and O^"*" - and therefore the total O abundance - from faint 
pure recombination lines. The results for NGC 2579 permit to extend our previous 
deter minations of the C, O and C/0 gas phase radial gradients of the inner Galactic 
disc (JEsteban et al.ll2005[ ) to larger galactocentric distances. We find that the chemi- 
cal composition of NGC 2579 is consistent with flatten gradients at its galactocentric 
distance. In addition, we have built a tailored chemical evolution model that repro- 
duces the observed radial abundance gradients of O, C and N and other observational 
constraints. We find that a levelling out of the star formation efficiency about and 
beyond the isophotal radius can explain the flattening of chemical gradients observed 
in the outer Galactic disc. 
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1 INTRODUCTION 

High-resolution spectroscopy with large aperture telescopes 
has permitted to measure very faint pure recombination 
lines (hereafter RLs) of heav y-element ions in G alactic and 
extragalactic H l l regions (seelGarcia-Roias fc Es tcban 2007; 



Peirnbertl l2003l : iLopez-Sanchez et al.1 l2007l : lEsteban et al.1 
20091 ). The brightest of these RLs are Cii 4267 A and 



* Based on observations collected at the European Southern Ob- 
servatory, Chile, proposal number ESO 382.C-0195(A). 
t E-mail: cel@iac.es 



those of multiple! 1 of O ll about 4650 A that can be used 
to determine the ionised gas phase C^""" and 0^+ abun- 
dances, respectively. In H ll regions, a systematic result is 
that abundances obtained with RLs are always between 1.5 
to 3 times larger (the so-called abundance discrepancy fac- 
tor, ADF) than those determined using the much brighter 
coUisionally excited lines (hereafter CELs) of the same ions. 
Cm] 1 909 A in the UV and those of [Cm] in the opt i- 
cal (e.g. lGarcia-Roias fc Estcban 2007: iEsteban et al.ll2009l ') . 
The origin of the abundance discrepancy (hereafter AD ) 
is still an unsolved problem. ITorres-Peimbert et al.l (|l980l ) 
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suggested that the presence of fluctuations in the spatial 
distribution of electron t emperature can produce suc h dis- 
crepancy but others as iTsamis fc Peguignotl (|2005l 'l and 
IStasiiiska et alj (|2007l ) consider that the AD may be pro- 
duced by small spatial scale chemical inho mogeneities in 
the interstellar medium. iTsamis et alJ (2011) have proposed 
that the presence of small and dense partially- ionised clumps 
along the line of sight axe affecting the CEL-derived abun- 
dances, producing unrealistic lower values with respect to 
the RL-b ased ones, and a simi l ar con clusion is found in the 
results bv lMesa-Delgado et alJ (|2012r ). More recently, a new 
hypothesis based on the presence of a K-distribution in the 
ener gy of the free electron s has been proposed to explain the 
AD l|Nichollset"aIll2012l ). 

Fro m deep echelle s pectr a of a sample of Galactic H ll 
regions, lEsteban et alJ (|2005l ) determined the carbon and 
oxygen gradients of the Galactic disc at galactocentric dis- 
tanc es between 6 . 3 and 10.4 kpc based on the intensity of 
RLs. ICarigi et al.l l|2005l ) built chemical evolution models in 
order to reproduce those abundance gradients as well as 
other observational constraints. They found the important 
result that only models including carbon yields that increase 
with metallicity for massive stars and decrease with metal- 
licity for low and intermediate mass stars can successfully 
reproduce simultaneously the gradients of both elements. 
Other result of that paper is that the fraction of carbon pro- 
duced by massive stars with respect to that produced by low 
and intermediate mass stars strongly depends on the age of 
the Galactic disc, as well as on the galactocentric distance. 
Therefore, it is clear that similar C and O abundance data of 
Hll regions at larger galactocentric distances are necessary 
to refine and extend chemical evolution models and explore 
possible changes across the Galactic disc. 

There are some evidences that gas phase abundance 
gradients ma y flatten out at the outer parts of the 
Galactic disc IfIcIi fc Silkevlll99ll:IVflchez fc Estebanll 19961: 



Maciel fc Quirezal Il999l : ICosta et all 12004 iMaciel et all 
2006|'l. However, this result is not supported by other stud- 



ies JDeharveng et al.ll2000l : lRudolph et al.ll2006l : lHenrv et all 
l20ld : ISalser et al.ll201lh . Chemical evolution models of the 
Galactic disc also obtain different solutions for the behaviour 
of the abund ance gradi ents at large galactocentric distances. 
For example. I Fu et al.l (|2009 ) find a st eepening of the gradi- 
ents a s radial distances increase and iMarcon-Uchida et al.l 
i|201CI ) can predict a fiattening or a steeping depending 
on the assumed behaviour of the star formation efficiency 
ac ross the Galactic disc. Moreover, chemodynamical models 
of ISamland et al.l (| 19971 ) also predict the development of a 
plateau in the outer regions of the Galaxy. Therefore, the 
finding of indisputable observational evidence of the pres- 
ence or not of such flattening is essential to increase our 
knowledge on important ingredients of chemical evolution 
models such as the timescale for disc formation and the den- 
sity threshold for star formation. 

NGC 2579 is a high surface brightness Galactic H ll re- 
gion that has been largely neglected because its misc las- 
sification as eithe r planetary nebula (I Acker et al.l 119921 ) or 
reflection nebula (|van den Bergh fc Herbstj 1973) and even 
its confusion with other ne arby objects JArchinal fc Hyned 
l2003l ). [Copetti et al] ((2007^ have brought into attention the 
interest of NGC 2579 presenting the flrst comprehensive ob- 
servational study of the nebular properties and stellar con- 




Figure 1. Combination of Ha and continuum image of 
NGC 2579. The position of the reference star DENIS J082054.8- 
361258 is indicated with a small circle. The grey (yellow) hori- 
zontal line centered at 5 arcsec to the north of the reference star 
indicates the area from which we extracted our one-dimensional 
spectra. North is up and east to the left. 



tent of this object. These authors derive consistent photo- 
metric and kinematic distances for NGC 2579 indicating 
that it is located at a galactocentric distance of 12.4 ± 0.7 
kpc (assuming the Sun at 8 kpc from the Galactic Centre). 
Therefore, NGC 2579 is probably the brightest Hll region 
located at that distant part of the Galaxy, well outside the 
solar circle. Another remarkable characteristic of this nebula 
is its relatively high ionization degree, which facilitates the 
determination of chemical abundances and the determina- 
tion of C'^"'' and O^^ abundances from RLs. Undoubtedly, 
NGC 2579 is an excellent probe for the exploration of the be- 
haviour of the ionised gas phase radial abundance gradients 
in the external parts of the Galactic disc. 

In ij2]we describe the observations and the data reduc- 
tion procedure. In i}3]we describe the emission line measure- 
ments and identifications as well as the reddening correction. 
In ^we present the physical conditions and ionic and total 
abundances determined for NGC 2579. In 3S]we describe the 
ingredients of our chemical evolution model and discuss the 
results in the light of the radial Galactic abundance gradi- 
ents. Finally, in iJB]we summarize our main conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

NGC 2579 was observed in service time on 2009 January 
23 at Cerro Paranal Observatory (Chile), using the UT2 
(Kueyen) of the Very Large Telescope (VL T) with the Ultra- 
violet Visual Echelle Spectrograph (UVES. IP'Odorico et al.l 
I2OO0I ). The standard settings of UVES were used covering 
the spectral range from 3570 to 10400 A. Some small spec- 
tral intervals could not be observed. These are: 5734-5833 
and 8425-8768 due to the physical separation between the 



© 2010 RAS, MNRAS 000, [THIS] 



NGC 2579 and the Galactic abundance gradients 3 




:"'""'"' 


1 1 1 . 


1 1 1 1 1 


7 on ON 


oil 


oil 7 


_ r 


n 




n _ 


*Ji* 


11 


m 


La 


1 1 1 1 1 



4865 



4370 



4640 4650 4660 4670 4660 
Wavelength (A) 



Figure 2. Sections of ttie UVES spectrum on NGC 2579 - un- 
corrected for reddening - siiowing the faint recombination lines 
of Cll 4267 A (left) and multiplet 1 of On about 4650 A (right). 



CCDs of the detector system of the red arm; and several 
much smaller wavelength ranges between 8800 to 10400 A 
because the last orders of the spectrum do not fit completely 
within the size of the CCD. 

Three consecutive exposures of 400 seconds - for the 
3570-3840 and 4820-6740 A ranges - and of 1800 seconds 
- for the 3770-4940 and 6735-10400 A ranges - each were 
added to obtain the final spectra. In addition, exposures of 
30 and 60 seconds were taken to obtain non-saturated flux 
measurements for the brightest emission lines. The full width 
at half-maximum (FWHM) at a given A was AA « A/8000. 
The slit was oriented east-west and the atmospheric disper- 
sion corrector (ADC) was used to keep the same observed 
region within the slit regardless of the air mass value. The 
slit width was set to 3" as a compromise between the spec- 
tral resolution needed and the desired signal-to-noise ratio 
of the spectra. The slit length was fixed to 10". The one- 
dimensional spectra we finally analysed were extracted for 
an area of 3" x 7'.'4. This area covers the brightest part of the 
nebula (see Figure[TJ, located 5" to the north of the star DE- 
NIS J082054.8-361258, a member of the stellar cluster that 
ionises the nebula. 

The spectra were reduced using the irae{j echelle reduc- 
tion package, following the standard procedure of bias sub- 
traction, aperture extraction, fiatfielding, wavelength cali- 
bration and flux calibrat ion. The standard star LTT 3218 
l|Hamuv et al.lll992l . 1 19941 ) was observed to perform the flux 
calibration. 



3 LINE INTENSITIES 

Line fluxes were measured by integrating all the fiux in the 
line between two given limits and over a local continuum es- 
timated by eye. In the case of line blending, a double Gaus- 
sian profile fit procedure was applied to measure the indi- 
vidual line intensities. All these measurements were made 
with the SPLOT routine of iraf. 

All line fiuxes of a given spectrum have been normalized 

^ IRAF is distributed by NOAO, which is operated by AURA, 
under cooperative agreement with NSF. 



to a particular bright emission line in each spectral range. 
For the bluest range (3570-3840 A), the reference line was 
H12 3750 A. For the 4820-6740 and 3770-4940 A ranges, H^ 
was used. For the reddest spectral interval (6735-10400 A), 
the reference was P20 8392 A. In order to produce a final 
homogeneous set of line flux ratios, all of them were rescaled 
to H/? considering the theoretical H12/H/? and P20/H/3 line 
ratios corresponding to the physical conditions of the gas. 

The spectral ranges present some overlapping at the 
edges. The final fiux of a line in the overlapping regions was 
the average of the values obtained in both adjacent spec- 
tra. A similar procedure was considered in the case of lines 
present in two consecutive spectral orders. The average of 
both measurements were considered as the adopted value of 
the line fiux. In all cases, the differences in the line flux mea- 
sured for the same line in different orders and/or spectral 
ranges do not show systematic trends and are always within 
the uncertainties. 

The identiflcation and laboratory wavelengths of the 
lines were obtained following our previous works on 
echelle spe ctroscopy of bright Galactic Hll regions (see 
iGarcia-Roias fc Esteban 2007 , and references therein). We 
also used a preliminar (unpublish ed) version of the spect ral 
synthesis code for nebulae X-SSN l|Pequignot et al.ll2012l ) to 
identify faint lines. This code was specially useful to identify 
lines in the reddest part of the spectrum, where confusion 
with telluric features is problematic. 

For a given line, the observed wavelength is determined 
by the centre of the baseline chosen for the flux integration 
procedures or the centroid of the line when a Gaussian fit is 
used (in the case of line blending). The final adopted values 
of the observed wavelength of a given line are relative to the 
heliocentric reference frame. 

The reddening coefficient, c(H/3), was determined from 
the comparison of the observed fiux ratio of the bright- 
est Balmer lines - Ha, H7, H(5 and Hz/ lines with re- 
spect to H/3 - and the th eoretical line ratios computed by 
IStorev fc Hummerl (|l995l ) for the physical conditions of the 
nebula. We have used t he reddening function , /(A), nor- 
malized to H;9 derived bv lCardelli et al.l l|l989t ) and assum- 
ing Rv = 3.1. The c(H/3) we obtain is 2.37±0.05, which 
is larger than the value s of th e reddening coefficient deter- 
mined bvlCopetti et al.l (|2007l ) who use the reddening func- 
tion bv lKalerl ( 19761 ) and obtain a range of values from 1.12 



to 1.65 depending on the position inside the nebula and the 
Balmer lines used. 

In Table [TJ the final list of line identifications (columns 
1-3) , observed wavelength corrected for heliocentric velocity 
(column 4), observed and dereddened fiux line ratios with 
respect to H/3 (columns 5 and 6) and the uncertainty of line 
ratios - in percentage - (last column) are presented. The 
quoted errors include the uncertainties in line flux measure- 
ment and error propagation in the reddening coefficient. 

In Figure (2] we show sections of our flux-calibrated 
echelle spectra showing the recombination lines of C 11 4267 
A and multiplet 1 of On around 4650 A. As it can be seen, 
these faint lines are well separated and show a remarkably 
high signal-to-noise ratio in our deep echelle spectrum. 
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Table 1. Observed, F(A), and reddening-corrected, /(A), line ratios [F{H/3) = 100] for NGC 2579. 



Ao 






^ohs 






Error 


(A) 


Ion 


Multiplet 


i^) 


F{X) 


7(A) 


i%r 


3587.28 


He I 


3 


3588.08 


0.039 


0.180 


17 


3613.64 


He I 


6 


3614.47 


0.067 


0.301 


12 


3634.25 


He I 


28 


3635.07 


0.052 


0.227 


14 


3662.26 


Hi 


H30 


3663.19 


0.045 


0.192 


16 


3663.40 


Hi 


H29 


3664.21 


0.036 


0.154 


18 


3664.68 


Hi 


H28 


3665.54 


0.056 


0.238 


14 


3666.10 


Hi 


H27 


3666.90 


0.059 


0.254 


13 


3667.68 


Hi 


H26 


3668.54 


0.076 


0.327 


12 


3669.47 


Hi 


H25 


3670.31 


0.080 


0.343 


11 


3671.48 


Hi 


H24 


3672.31 


0.088 


0.377 


11 


3673.76 


Hi 


H23 


3674.60 


0.090 


0.384 


11 


3676.37 


Hi 


H22 


3677.20 


0.114 


0.483 


10 


3679.36 


Hi 


H21 


3680.21 


0.156 


0.663 


8 


3682.81 


Hi 


H20 


3683.62 


0.138 


0.581 


9 


3686.83 


Hi 


H19 


3687.68 


0.126 


0.531 


9 


3691.56 


Hi 


HIS 


3692.42 


0.186 


0.779 


8 


3697.15 


Hi 


H17 


3697.98 


0.250 


1.042 


7 


3703.86 


Hi 


H16 


3704.71 


0.262 


1.085 


7 


3705.04 


He I 


25 


3705.87 


0.121 


0.501 


9 


3711.97 


Hi 


H15 


3712.81 


0.325 


1.338 


6 


3721.83 


[Siii] 


2F 


3722.68 


0.646 


2.631 


6 


3721.94 


Hi 


H14 










3726.03 


[On] 


IF 


3726.90 


16.00 


64.94 


5 


3728.82 


[On] 


IF 


3729.65 


12.51 


50.63 


5 


3734.37 


Hi 


H13 


3735.22 


0.467 


1.880 


6 


3750.15 


Hi 


H12 


3751.00 


0.690 


2.737 


6 


3770.63 


Hi 


Hll 


3771.49 


0.815 


3.168 


5 


3797.90 


Hi 


HIO 


3798.76 


1.179 


4.462 


5 


3819.61 


He I 


22 


3820.49 


0.262 


0.969 


7 


3833.57 


He I 


62 


3834.47 


0.017 


0.062 


28 


3835.39 


Hi 


H9 


3836.26 


1.710 


6.229 


5 


3856.02 


Sin 


IF 


3856.90 


0.026 


0.093 


22 


3856.13 


On 


12 










3862.59 


Sin 


1 


3863.49 


0.022 


0.079 


24 


3867.49 


He I 


20 


3868.34 


0.042 


0.150 


16 


3868.75 


[Ne III] 


IF 


3869.64 


3.409 


11.99 


5 


3871.82 


He I 


60 


3872.68 


0.018 


0.062 


27 


3888.65 


He I 


2 


3889.85 


3.911 


13.46 


5 


3889.05 


Hi 


H8 










3918.98 


Cii 


4 


3919.82 


0.008 


0.028 




3920.68 


On 


4 


3921.50 


0.023 


0.076 


23 


3926.53 


He I 


58 


3927.42 


0.035 


0.114 


18 


3964.73 


He I 


5 


3965.63 


0.289 


0.914 


6 


3967.46 


[Ne III] 


IF 


3968.37 


1.131 


3.567 


5 


3970.07 


Hi 


H7 


3970.97 


4.912 


15.45 


5 


4009.22 


He I 


55 


4010.19 


0.067 


0.203 


12 


4026.21 


He I 


18 


4027.12 


0.733 


2.161 


5 


4068.60 


[Sii] 


IF 


4069.56 


0.239 


0.670 


7 


4069.62 


On 


10 


4070.66 


0.033 


0.092 


19 


4072.15 


On 


10 


4073.05 


0.022 


0.063 


24 


4075.86 


On 


10 


4076.81 


0.021 


0.058 


25 


4076.35 


[Sii] 


IF 


4077.32 


0.075 


0.209 


12 


4101.74 


Hi 


H6 


4102.67 


9.500 


25.60 


5 


4120.82 


He I 


16 


4121.72 


0.085 


0.225 


11 


4143.76 


He I 


53 


4144.72 


0.136 


0.347 


8 


4153.3 


On 


19 


4154.29 


0.012 


0.030 


36 


4168.97 


He I 


52 


4169.96 


0.019 


0.048 


26 


4267.15 


Cii 


6 


4268.14 


0.071 


0.157 


12 


4287.40 


[Fen] 


7F 


4288.40 


0.013 


0.027 


34 


4317.14 


On 


2 


4318.07 


0.010 


0.021 


39 


4326.4 


Oi 




4327.45 


0.010 


0.021 


38 
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Table 1. continued 



Ao 






Aofes 






Error 


(A) 


Ion 


Multiple! 


(A) 


F{\) 


/(A) 


(%)- 


4340.47 


Hi 


H7 


4341.45 


23.67 


47.27 


4 


4345.56 


On 


2 


4346.57 


0.019 


0.038 


26 


4349.43 


On 


2 


4350.41 


0.020 


0.039 


25 


4361.60 


Nil 




4362.55 


0.012 


0.023 


35 


4363.21 


[Oiii] 


2F 


4364.19 


0.968 


1.877 


5 


4366.89 


On 


2 


4367.90 


0.021 


0.041 


24 


4368.19 


Oi 


5 


4369.32 


0.019 


0.037 


26 


4368.25 


Oi 


5 










4387.93 


He I 


51 


4388.93 


0.306 


0.574 


6 


4437.55 


He I 


50 


4438.56 


0.043 


0.076 


16 


4471.48 


He I 


14 


4472.52 


2.727 


4.592 


4 


4607.13 


[Fein] 


3F 


4608.13 


0.023 


0.032 


23 


4609.44 


On 


93 


4610.42 


0.009 


0.012 




4630.54 


Nil 


5 


4631.59 


0.026 


0.035 


22 


4638.86 


On 


1 


4639.90 


0.037 


0.051 


17 


4641.81 


On 


1 


4642.86 


0.054 


0.073 


14 


4649.13 


On 


1 


4650.22 


0.071 


0.094 


12 


4650.84 


On 


1 


4651.87 


0.034 


0.045 


18 


4658.10 


[Fe III] 


3F 


4659.18 


0.417 


0.547 


5 


4661.63 


On 


1 


4662.64 


0.044 


0.057 


16 


4667.01 


[Felll] 


3F 


4668.04 


0.017 


0.022 


28 


4673.73 


On 


1 


4674.79 


0.009 


0.012 




4676.24 


On 


1 


4677.34 


0.025 


0.032 


22 


4701.53 


[Felll] 


3F 


4702.65 


0.116 


0.143 


9 


4711.37 


[Ariv] 


IF 


4712.49 


0.018 


0.023 


27 


4713.14 


He I 


12 


4714.25 


0.501 


0.611 


5 


4733.93 


[Fe III] 


3F 


4734.99 


0.042 


0.050 


16 


4740.16 


[Ariv] 


IF 


4741.32 


0.016 


0.019 


29 


4754.83 


[Felll] 


3F 


4755.82 


0.083 


0.096 


11 


4769.60 


[Fe III] 


3F 


4770.55 


0.044 


0.050 


15 


4777.88 


[Felll] 


3F 


4778.81 


0.021 


0.023 


25 


4861.33 


Hi 


H/3 


4862.42 


100.0 


100.0 


4 


4881 


[Felll] 


2F 


4882.15 


0.180 


0.176 


7 


4889.62 


[Fell] 


4F 


4890.82 


0.014 


0.014 


31 


4894.59 


Oi 




4895.86 


0.010 


0.010 


38 


4921.93 


He I 


48 


4923.05 


1.202 


1.109 


4 


4924.5 


[Fciii] 


2F 


4925.67 


0.024 


0.022 


22 


4931.32 


[Oiii] 


IF 


4932.36 


0.045 


0.041 


15 


4958.91 


[Oiii] 


IF 


4960.09 


139.9 


122.8 


4 


4985.9 


[Fe III] 


2F 


4986.97 


0.072 


0.061 


12 


4987.2 


[Fe III] 


2F 


4988.38 


0.079 


0.067 


11 


5006.84 


[Oiii] 


IF 


5008.03 


432.8 


356.2 


4 


5011.30 


[Felll] 


IF 


5012.46 


0.063 


0.052 


12 


5015.68 


He I 


4 


5016.86 


2.802 


2.279 


4 


5041.03 


Sin 


5 


5042.21 


0.100 


0.079 


10 


5055.98 


Sin 


5 


5057.21 


0.172 


0.133 


7 


5191.82 


[Arm] 


3F 


5192.91 


0.100 


0.064 


10 


5197.90 


[Ni] 


IF 


5199.27 


0.149 


0.095 


8 


5200.26 


[Ni] 


IF 


5201.69 


0.092 


0.059 


10 


5261.62 


[Fell] 


19F 


5263.00 


0.022 


0.013 




5270.40 


[Felll] 


IF 


5271.75 


0.341 


0.198 


6 


5412.00 


[Fe III] 


IF 


5413.40 


0.038 


0.018 


17 


5517.71 


[CI III] 


IF 


5518.98 


0.961 


0.415 


5 


5537.88 


[CI III] 


IF 


5539.13 


0.914 


0.387 


5 


5875.64 


He I 


11 


5877.03 


32.51 


14.44 


4 


5941.65 


Nil 


28 


5943.18 


0.036 


0.015 


17 


5978.93 


Sin 


4 


5980.31 


0.141 


0.057 


8 


6300.30 


[Oi] 


IF 


6301.91 


1.164 


0.359 


4 


6312.10 


[Sni] 


3F 


6313.55 


5.293 


1.615 


4 


6347.11 


Sin 


4 


6348.58 


0.505 


0.150 


5 


6363.78 


[Oi] 


IF 


6365.40 


0.444 


0.130 


5 


6371.36 


Sin 


2 


6372.87 


0.260 


0.076 


6 
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Table 1. continued 



Ao 






^obs 






Error 


(A) 


Ion 


Multiplet 


(A) 


F{\) 


7(A) 


{%)" 


6548.03 


[Nil] 


IF 


6549.61 


30.40 


7.681 


4 


6562.82 


Hi 


Ha 


6564.31 


1161 


290.0 


4 


6578.05 


Cii 


2 


6579.60 


0.599 


0.148 


5 


6583.41 


[Nil] 


IF 


6585.02 


96.78 


23.79 


4 


6678.15 


He I 


46 


6679.77 


15.34 


3.508 


4 


6716.47 


[Sii] 


2F 


6718.03 


15.70 


3.487 


4 


6730.85 


[Sii] 


2F 


6732.60 


18.66 


4.101 


4 


6734.08 


He I 


1/22 


6735.84 


0.025 


0.005 


22 


6734.42 


[Feiv] 












6739.8 


[Feiv] 




6741.65 


0.014 


0.004 


31 


6744.1 


He I 


1/22 


6745.93 


0.011 


0.003 


37 


6755.95 


He I 


1/20 


6757.60 


0.013 


0.004 


33 


6755.99 


[Feiv] 












6785.68 


He I 


1/18 


6787.67 


0.020 


0.005 


25 


6787.21 


Cii 


14 


6788.84 


0.010 


0.003 




6791.47 


Cii 


14 


6793.15 


0.011 


0.003 


37 


6804.95 


He I 


1/17 


6806.75 


0.025 


0.007 


22 


6855.88 


He I 


1/15 


6857.73 


0.023 


0.006 


23 


6933.99 


He I 


1/13 


6935.78 


0.031 


0.008 


19 


6989.53 


He I 


1/12 


6991.31 


0.051 


0.012 


14 


7062.34 


He I 


1/11 


7064.10 


0.085 


0.019 


11 


7065.28 


He I 


10 


7066.90 


26.66 


5.885 


4 


7135.78 


[Ar III] 


IF 


7137.45 


65.90 


13.87 


4 


7155.16 


[Fell] 


14F 


7157.11 


0.088 


0.018 


10 


7160.56 


He I 


1/10 


7162.42 


0.135 


0.028 


8 


7231.32 


Co 


3 


7233.25 


0.180 


0.036 


7 


7236.42 


Cii 


3 


7238.15 


0.508 


0.100 


5 


7281.35 


Ho I 


45 


7283.06 


3.512 


0.671 


4 


7298.04 


He I 


1/9 


7299.84 


0.102 


0.019 


10 


7318.39 


[On] 


2F 


7320.77 


4.099 


0.765 


4 


7319.99 


[On] 


2F 


7321.85 


13.18 


2.455 


4 


7329.66 


[On] 


2F 


7331.41 


6.969 


1.291 


4 


7330.73 


[On] 


2F 


7332.49 


7.297 


1.351 


4 


7377.83 


[Niii] 


2F 


7379.90 


0.123 


0.022 


9 


7423.64 


Ni 


3 


7425.70 


0.022 


0.004 


24 


7442.3 


Ni 


3 


7444.39 


0.070 


0.012 


12 


7452.54 


[Fell] 


14F 


7454.55 


0.033 


0.006 


19 


7499.85 


He I 


1/8 


7501.77 


0.238 


0.040 


6 


7504.94 


On 




7506.86 


0.020 


0.0033 


25 


7513.33 


He I 




7515.25 


0.017 


0.0028 


28 


7519.29 


He I 




7521.52 


0.036 


0.006 


17 


7519.49 


Cii 


16.8 










7519.93 


Cn 


16.8 










7530.54 


[01 IV] 


IF 


7532.51 


0.027 


0.0043 


21 


7530.57 


Oil 


16.8 










7751.10 


[Ar III] 


2F 


7752.92 


20.94 


3.001 


4 


7774.17 


Oi 


1 


7776.17 


0.031 


0.0044 


19 


7775.39 


Oi 


1 


7777.32 


0.017 


0.0024 


28 


7816.13 


He I 


1/7 


7817.98 


0.402 


0.055 


4 


8015.96 


He I 


4/20 


8017.94 


0.024 


0.0029 


23 


8035.06 


He I 


4/19 


8037.14 


0.038 


0.0047 


17 


8045.62 


[01 IV] 


IF 


8047.82 


0.59 


0.0072 


13 


8057.55 


He I 


4/18 


8059.50 


0.042 


0.0051 


16 


8094.11 


He I 


4/10 


8096.18 


0.039 


0.0046 


17 


8116.44 


He I 


4/16 


8118.45 


0.049 


0.0058 


14 


8203.86 


Ho I 


4/14 


8205.98 


0.113 


0.0126 


9 


8210.72 


Ni 


2 


8213.11 


0.045 


0.0051 


15 


8216.34 


Ni 


2 


8218.82 


0.078 


0.0087 


11 


8223.14 


Ni 


2 


8225.21 


0.206 


0.023 


7 


8233.21 


Hi 


P50 


8235.13 


0.090 


0.0100 


10 


8234.43 


Hi 


P49 


8236.39 


0.114 


0.0125 


9 


8235.74 


Hi 


P48 


8237.71 


0.113 


0.0124 


9 
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Table 1. continued 



Ao 






^ohs 






Error 


(A) 


Ion 


Multiple! 


(A) 


F{X) 


/(A) 


{%r 


8237.13 


Hi 


P47 


8239.05 


0.183 


0.020 


7 


8238.61 


Hi 


P46 


8240.55 


0.202 


0.022 


7 


8240.19 


Hi 


P45 


8242.02 


0.283 


0.031 


6 


8241.88 


Hi 


P44 


8244.12 


0.304 


0.033 


6 


8243.69 


Hi 


P43 


8245.60 


0.274 


0.030 


6 


8245.64 


Hi 


P42 


8247.56 


0.379 


0.042 


5 


8247.73 


Hi 


P41 


8249.66 


0.346 


0.038 


5 


8249.2 


Hi 


P40 


8251.89 


0.343 


0.038 


5 


8252.4 


Hi 


P39 


8254.35 


0.404 


0.044 


5 


8255.02 


Hi 


P38 


8257.04 


0.319 


0.035 


6 


8257.85 


Hi 


P37 


8259.82 


0.367 


0.040 


5 


8260.93 


Hi 


P36 


8262.84 


0.460 


0.050 


5 


8264.28 


Hi 


P35 


8266.28 


0.660 


0.072 


5 


8265.7 


He I 


2/9 


8267.71 


0.055 


0.0059 


14 


8267.94 


Hi 


P34 


8269.89 


0.555 


0.060 


5 


8271.93 


Hi 


P33 


8273.80 


0.512 


0.055 


5 


8276.31 


Hi 


P32 


8278.33 


1.543 


0.167 


4 


8281.12 


Hi 


P31 


8283.03 


0.865 


0.093 


5 


8290.56 


He I 


6/27 


8292.71 


0.025 


0.0027 


22 


8292.31 


Hi 


P29 


8294.28 


0.597 


0.064 


5 


8298.83 


Hi 


P28 


8300.82 


0.785 


0.084 


5 


8306.11 


Hi 


P27 


8308.06 


1.067 


0.114 


4 


8314.26 


Hi 


P26 


8316.22 


1.055 


0.112 


4 


8323.42 


Hi 


P25 


8325.39 


1.290 


0.136 


4 


8329.86 


He I 


6/23 


8331.76 


0.038 


0.0040 


17 


8333.78 


Hi 


P24 


8335.76 


1.396 


0.146 


4 


8345.55 


Hi 


P23 


8347.51 


1.599 


0.167 


4 


8359.00 


Hi 


P22 


8360.93 


2.018 


0.209 


4 


8361.73 


He I 


1/6 


8363.81 


0.822 


0.085 


5 


8374.48 


Hi 


P21 


8376.43 


1.994 


0.205 


4 


8376.56 


He I 


6/20 


8378.62 


0.066 


0.0068 


12 


8392.40 


Hi 


P20 


8394.36 


2.325 


0.237 


4 


8413.32 


Hi 


P19 


8415.24 


4.827 


0.487 


4 


8421.97 


He I 


6/18 


8424.04 


0.118 


0.012 


9 


8862.79 


Hi 


Pll 


8862.50 


11.71 


1.410 


4 


8914.77 


He I 


2/7 


8917.12 


0.119 


0.014 


9 


8930.77 


He I 


10/11 


8933.30 


0.044 


0.0051 


15 


8996.99 


He I 


6/10 


8998.67 


0.489 


0.056 


5 


9014.91 


Hi 


PIO 


9017.24 


7.810 


0.882 


4 


9063.32 


He I 


4/8 


9065.61 


0.466 


0.052 


5 


9068.90 


[Siii] 


IF 


9071.05 


198.7 


21.94 


4 


9123.6 


[CI 11] 


IF 


9126.03 


0.131 


0.014 


8 


9210.28 


He I 


6/9 


9212.52 


0.635 


0.066 


5 


9213.24 


He I 


7/9 


9215.60 


0.239 


0.025 


6 


9229.01 


Hi 


P9 


9231.26 


19.40 


2.008 


4 


9262.67 


Oi 


8 


9265.07 


0.021 


0.0021 


25 


9265.94 


Oi 


8 


9268.33 


0.031 


0.0031 


19 


9463.57 


He I 


1/5 


9465.78 


1.120 


0.106 


4 


9526.17 


He I 


6/8 


9528.41 


0.465 


0.043 


5 


9530.60 


[Siii] 


IF 


9532.95 


520.5 


48.02 


4 


9545.97 


Hi 


P8 


9548.06 


27.04 


2.48 


4 


9625.70 


He I 


6/8 


9628.23 


0.207 


0.019 


7 


9682.39 


He I 


9/8 


9685.01 


0.041 


0.0036 


16 


9701.87 


[Feiii] 




9704.40 


0.131 


0.0113 


8 


9702.65 


He: 


75 


9705.29 


0.172 


0.015 


7 


9903.46 


Cii 


17.2 


9906.40 


0.389 


0.032 


6 


10027.70 


He I 


6/7 


10030.66 


2.050 


0.159 


4 


10031.2 


He I 


7/7 


10033.94 


0.657 


0.051 


5 


10049.40 


Hi 


P7 


10051.76 


57.70 


4.450 


4 


10138.42 


He I 


10/7 


10140.91 


0.247 


0.019 


6 



Error of the dereddened flux ratios. Colons indicate errors larger than 40 per cent. 
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4 PHYSICAL CONDITIONS AND CHEMICAL 
ABUNDANCES OF NGC 2579 

The physical conditions of the ionised gas: electron tem- 
perature, Tc, and density, rio, have been derived from 
the usual CEL rat ios, usins the iraf task temden of the 
package nebular (IShaw fc Dufour 1995) with the atomic 
dataset compiled by iGarcia-Roias et alj l|2005t) . Electron 
densities have been derived from [N l] AA5198/5200, [O ll] 
AA3726/3729, [Sii] AA6717/6731, [Cliii] AA5518/5538, and 
[Ariv] AA4711/4740 hue ratios. We have also determined 
rie from the line ratios of several [Fe ill] lines fo llowing the 
procedure described bv lGarcia-Roias et alJ (|2006l ). Electron 
temperatures have been calculated using several sets of auro- 
ral to nebular line intensity ratios: [O ll] AA7319, 7330/ [O ll] 
AA3726, 3729, [S ii] AA4069, 4076/[Sii] AA6717, 6731, [Qui] 
A4363/[Oiii] AA4959, 5007, [S iii] A6312/[Siii] AA9069, 9532 
and [Arm] A5192/[Ar iii] AA7136, 7751. Unfortunately, it 
was not possible to calculate To ([Nil]) because the auroral 
[N ll] A5755 A line was in an incomplete spectral order. 

The procedure for the determination of the physical 
conditions was the following: an initial To-value of 8,000 
K was assumed in order to derive a first approximation to 
the different Uc determinations; then these preliminary n^- 
values were used to recompute Tc, and finally, we iterated 
until convergence to obtain the final adopted values of n^ 
and Tc. ro([Oll]) has been corrected from the contribution 
to [O ll] AA7319, 73 30 due to recom bination following the 
formulae derived bv lLiu et al.l (|2000f ). This correction only 
contributes a 2.4% to the intensity of [On] AA7319, 7330 
lines. The physical conditions determined for NGC 2579 are 
shown in Table [21 

NGC 257 9 pres e nt a n on-uinform density structure. 
ICopetti et all (|2000l . 120071 ') showed that the no([Sii]) 
presents a strong spatial variation, with the density rang- 
ing from about 1800 cm~"^ at the brightest eastern-central 
areas to less than 100 cm~^ at the outer parts of the neb- 
ula. This density variation was interpreted as steep ra- 
dial gradient, similar to that found in the Orion Neb- 
ula (|Osterbrock fc Flatherlll959l ). However, the density esti- 
mates listed in Tableware rather sim ilar and between 1000 
and 2000 cm~^. ICopetti et al.l l|2007h obtained nc([Sll]) ^ 
1200 cm~^ at the position of the present observations and 
an upper limit of ne([Cllll]) of about 2100 cm~^ from the in- 
tegrated spectrum. Those authors derive values of re([Nll]) 
and Te([Olll]) of 11000 and 9000 K, respectively, consistent 
with our determinations. 

Ionic abundances of N+, 0+, 0^+, ^e^+ , S+, S^+, C1+, 
Cl^+ , Cl^+ , Ar2+ and Ar^+ have been derived from CELs 
under the two-zone scheme and t^ — 0, using the nebu- 
lar package, except in the case of Cl^ for which we have 
used another analysis package (see below) . We have assumed 
Te([Oll]) for the abundance calculations of low ionization 
potential ions: N+, O"*" and S"*"; and rc([Olll]) for the high 
ionization potential ones: O^"*", Ne^^, S^"^, Cl'^'*', Ar'^"'" and 
Ar^+. We have adopted rie = 1360 ± 170 for all the ions. This 
value is the weighted mean of the densities obtained from 
[S ii] , [O ii] and [CI iii] line rati os. For Ne^^ we have used the 
updated atomic data listed in lCarcia-Roias et al.l ([20091) in- 
stead of those included by default in nebular. The older 
atomic dataset gave an inconsistency of about 0.15 dex in 
the abundances obtained from the two [Ne ill] lines observed. 



Table 2. Physical Conditions of NGC 2579. 





Lines 


Value 


nc (cm ^) 


[Ni] 


1500±700 




[On] 


1570±280 




[Sii] 


1140±230 




[Feni] 


2330±800 




[CI III] 


1900±600 




[Ariv] 


1670if~ 


T. (K) 


[On] 


10450±250 




[Sii] 


7840±360 




[Oiii] 


9410±160 




[Siii] 


10760±230 




[Arm] 


8650±280 



Many [Felll] lines have been detected in the spectrum of 
NGC 2579. We have calculated the Fe^+/H+ ratio using 
a 34- level model a tom that includes the collision strengths 
from lZhand (|l996r ). transition prob abilities of|Quineti fil996l ) 
as well as the transitions found bv lJohansson et al.l ((20001). 
The average value of Fe^^ abundance has been obtained 
from 12 individual emission lines and assuming Tc([Oll]) 
as the representative temperature for this ion. We have de- 
tected [Fell] 4287 A, the only line of this ion present in 
our spectrum. The intensity of this li ne is affected by con- 
tinuum or starlight fiuorescence (see 'Rodriguez 19991) and 
therefore we can not determine an accurate abundance from 
this line. Unfortunately, the brightest and less sensitive to 
fiuorescence [Fell] line, [Fell] 8616 A is in one of our ob- 
servational gaps. Several [Feiv] lines have been detected in 
NGC 2579 however only [Feiv] 6739.8 A is not blended with 
other lines. Using that single line, the Fe^^/H^ ratio has 
been derived using a 33-level model at om where all col- 
lision strengths are those calculated by IZhang fc PradhanI 
(Il997h and the transition probabilities recommended by 
iFroese Fischer et al.l (|2008l '). We have asummed Tc([Olll]) 
to derive the Fe^"*" abundances. The Cl^ abundance cannot 
be deri ved from the nebular routines, instead we have used 
PYNEB (jLuridiana et al.ll2012l ). which is an updated version 
of the NEBULAR package written in python programming 
language. For CI'*' we h ave used the line wavelengths an d en- 
ergy levels obtained bv lRadziemski fc KaufmanI (119741 ) . the 



transition probabilities compiled bv | Mendozal ( 19831 ) and the 
collision strengths from iTavall (|200J) . Ionic abundances are 
presented in Table |31 

We have measured several He I emission lines in the 
spectrum of NGC 2579. These lines arise basically from re- 
combination but they can be affected by collisional exci- 
tation and self-absorption effect s. We have used the effec- 
tive recombination coefficients of IPorter et al . (2005) - with 
the interpolation formulae pr ovided by Porter et al.l (120071) 
- for He I as well as those of IStorev fc Hummed ( 19951 ) for 
Hi, in order to calculate the He"*" a bundance. The coUi- 
sional contribution was estimated from lSawev fc BerringtonI 
(1993) and Kingdon & Fcrland (1995|), and the optical 
depth in the triplet lines were derived from the computa- 
tions by iBeniamin et al.l (|2002l ). We have determined the 
He"'"/H"'" ratio from a max imum likelihood method (MLM, 
iPeimbert eraLll2000l . [200^ . 

To determine the He'''/H+ ratio, re(Hel), the tem- 
perature fiuctuations parameter (|Peimbertl 119671 ) for He I, 
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Table 3. Ionic and total abundances^ 





This work 


C07b 




i2 = 0.000 


t2 = 0.045 
± 0.007 


i2 = 0.000 


Ionic abundances from CELs*^ 


N+ 


6.62±0.02 


6.75±0.03 


6.63±0.12 


0+ 


7.70±0.03 


7.84±0.03 


7.65±0.16 


02+ 


8.19±0.02 


8.46±0.03 


8.31±0.06 


Ne2+ 


7.23±0.03 


7.53±0.06 


7.59±0.08 


S+ 


5.24±0.02 


5.36±0.03 


5.30±0.12 


g2+ 


6.54±0.02 


6.84±0.06 


6.78±0.08 


C1+ 


3.38±0.04 


3.49±0.06 


- 


C12+ 


4.80±0.02 


5.06±0.05 


5.03±0.05 


Cl3+ 


2.78±0.07 


3.01±0.10 


— 


Ar2+ 


6.14±0.02 


6.37±0.04 


- 


Ar3+ 


3.66±0.09 


3.94±0.10 


- 


Fe2+ 


5.27±0.10 


5.40±0.10 


- 


Fe3+ 


5.14±0.29 


5.37±0.30 


- 


Ionic abundances from RLs'* 


He+ 


10.94±0.01 


10.96±0.04 


C2+ 


8.13±0.05 


— 


0+ 


7.85±0.07 


— 


02+ 


8.46±0.03 


- 


Total abundances from CELs'^ 


N 


7.31±0.04 


7.54±0.05 


7.37±0.06 


O 


8.31±0.02 


8.55±0.03 


8.39±0.05 


Ne 


7.35±0.02 


7.62±0.07 


7.68±0.06 


S 


6.64±0.02 


6.93±0.06 


6.91±0.11 


CI 


4.82±0.02 


5.08±0.05 


5.16±0.08 


Ar 


6.17±0.02 


6.38±0.04 


- 


Fe 


5.74±0.08 


5.94±0.08 


- 


Total abundances from RLs'^ 


He 


10.96±0.01 


10.96±0.04 


C 


8.30±0.05 


— 


O 


8.56±0.03 


- 



" In units of 12+log(X+n/H+) . 
' Copetti et al. (2007). 
'^ CELs: CoUisionally excited lines. 
'^ RLs: Recombination lines. 



t'^(Hel), and the optical depth in the He I A3889 line, rsggg, 
in a self-consistent manner, we have used the adopted den- 
sity obtained from the CEL ratios (see Table [2|), and a set 
of 16 /(Hel)//(Hl) line ratios (AA 3614, 3819, 3889, 3965, 
4026, 4121, 4388, 4471, 4713, 4922, 5016, 5048, 5876, 6678, 
7065, 7281 A). We have obtained the best sets of values of 
the three unknowns by minimizing x^- The lowest x^ pa- 
rameters we obtain are in the range from 14 to 28, which 
indicate reasonably good fits. The final adopted He"*" abun- 
dance is 12-Hlog(He+/H+) = 10.94 ± 0.01 (see Table EJ. 
Interestingly, the range of values of f^(Hel) we obtain in 
these fits is between 0.025 and 0.065, consistent with that 
obtained assuming that the ADF(0^''') is produced by the 
presence of temperature fluctuations in the ionised gas (see 
below) . 

We have detected several Cll lines in our spectrum, 
but most of them are produced by re sonance fluorescence 
by starlight (JEsteban et al.ll 19981 . [200I ). Only Cll 4267 and 



9903 A are pure RLs and permit to derive the C^^/H+ ra- 
tio. Two RLs of multiplet 1 of Ol are detected in the red 
part of the spectrum, namely O I 7774 and 7775 A, and for- 
tunately these important lines are not affected by telluric 
emission. We also detect and measure a large number of 
RLs of On. In particular, seven lines of multiplet 1, four of 
multiplet 2 and three of multiplet 10 that are produced by 
pure recombination and can be used to derive the O^^/H"*" 
ratio. We have use d re([Oin]) and the effective reco mbina- 
tion c oefficients of iDavev et al.l (l2000|) for C 11 and IStorevI 
(|l994l ) for 11 to determine the C^+ and 0^+ abundances; 
an d Te([On]) and the effect ive r ecombination coefficients 
of lEscalante fc VictoJ j 19921 ) and iPeeauignot et al.l (|l99ll ) 
for the O"*" abundance. The C^'^ abundance obtained from 
the Cll 4267 and 9903 A lines are 12 + log (C^+/H+) = 
8.18 ± 0.05 and 8.10 ± 0.03, respectively We adopted the 
weighted mean of both values for the final C^^ abundance 
included in Table O In the case of the O^/H"'" we adopted 
the mean value of the abundances obtained from each of 
the two lines detected and the two sets of effective recom- 
bination coefficients indicated above. Following our usual 
methodology to minimize uncertainties, we have derived the 
O^"*" ab undance from the est imated total flux of each multi- 
plet fsee lEsteban et al.lll998l ). In particular, we obtain values 
of 12-Mog(0^+/H+) = 8.46±0.03, 8.51±0.06 and 8.40±0.07 
for multiplets 1, 2 and 10 respectively, remarkably similar 
and consistent within the errors. We have assumed the abun- 
dance obtained from multiplet 1 as the representative one for 
O^^ (a weighted mean value of the abundances given by each 
multiplet gives almost the same value). Our spectrum also 
shows a relatively large number of permitted lines of other 
heavy-element ions (Nil, Ol and Sill) but they are not pro- 
duced by pure recombinatio n and we can not derive re liable 
abundances from them fsee lEsteban et al.lll998l . |2004| ). 

Our deep spectrum permits to calculate O^ and O^^ 
abundances from two kinds of lines - RLs and CELs. As for 
all H 11 regions where both kinds of lines ha ve been observe d 
(e.g. Garcia-Roias & Estcban 2007; Estc ban et all |2009| ). 
the abundance determined from RLs is larger than that 
obtained from CELs. Then, we can compute the so-called 
abundance discrepancy factor for O"*" and O^^, ADF(X') , 
which is deflned in its logarithmic form as: 



ADF(X') = log(X7H+)RLs - log(X'/H+)cELs. 



(1) 



The value of the ADF(0^+) we obtain for NGC 2579 is 0.27 
± 0.03, remarkably similar to the values of this quantity 
found for other Galactic an d extragalactic H 11 re gions, which 
mean value is 0.26 ± 0.09 (jEsteban et al.ll2009l ). In the case 
of 0+, we obtain ADF(0+) = 0.15 ± 0.08, which is also in 
excellent agreement with the values of between 0.15 and 0.20 
we obtain for other well observed Galactic Hll regions as 
M 8, M 20 and the Orion nebula IjGarcia-Roias fc EstebanI 
l2007h . 

If - as an hypothesis - we assume the validity 
of the temperature fluctuations paradigm and that this 
phenomenon produces the a bundance discrepancy (see 
iGarcia-Roias fc EstebanI 120071 ). we can estimate the t^ pa- 
rameter that produces the agreement between the abun- 
dance of O^"*" determined from CELs and RLs, which re- 
suhs to be t^ = 0.045 ± 0.007. In Table Owe include ionic 
abundances determined for f^ = - the standard procedure 
considering no temperature fluctuations - and assuming t^ 
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= 0.045 ± 0.007. These last calculati ons have been made 
follow ing the formalism outlined by IPeimbert fc Costerd 
1 19691 ). It is remarkable that the t^ parameter value we ob- 
tain from the comparison of O^ abundances determined 
from CELs and RL is 0.048 ± 0.029, in excellent agree- 
ment with the value we obtain for O^^. This consistency 
supports that the phenomenon that produce the abundance 
discrepancy and temperature fluctuations may be connected 
or even be the same, at least in Hll regions (see discussion 
in lGarcia-Roias fc Estebanll2007^ . 

We have adopted a set of ionization correction factors 
(ICFs) to correct for the unseen ionization stages and de- 
rive the total gas phase abundances of the different ele- 
ments, except in the cases of O and CI, for which we have 
measured emission lines of all the expected ionic species of 
these elements. The final total abundances for t^ = and 
t^ = 0.045 ± 0.007 are presented in Table [3] For He, C, 
N, S and Ne we have adopted the same ICF schemes used 
bv lGaJcia-Roias fc EstebanI (|2007l ) in order to facilitate the 
comparison with data of other Galactic Hll regions and 
study the abundance gradients. The total helium abundance 
has been corrected for the presence of neut r al hel ium using 
the expression proposed bv lPeimbert et al.l l| 19921 ) based on 
the similarity of the ionization potentials of He° and S""" . In 
the case of C we have ad opted the ICF f C"*") derived from 
photoionization models of iGarnett et al.l l| 19991 ). In order to 
derive the total abundance of nitrogen we have use d the 
ICF derived from the models bv lMathis fc Rosal (|l99lh : note 
that these ICF values yield N abundances about 0.11 dex 
higher than those obtained from the usual formulation by 
IPeimbert fc Costerd (|l969l ) based on the similarity of the 
ionization potential of N^ and O"'' . The total abundance of 
oxygen is calculated as the sum of O"'" and O^''' abundances. 
The only measurable CELs of Ne in the optical range are 
those of Ne^+ but the fraction of Ne""" may not be negli- 
gib le in the nebula. We have adopted the usual expression 
of IjPeimbert fc Costerdll969l ) to obtain the total Ne abun- 
dance. This scheme seems to be apropriate for the ioniza- 
tion degree of NGC 2579. We have measured CELs of two 
ionizatioi i stages of S: S^ and S^^, and used the ICF pro- 
posed bv lStasiriskal l| 19781 ) to take into account the presence 
of some S^"^. Chlorine shows lines of the three ionization 
stages we expect in NGC 2579. Its total abundance is sim- 
ply the sum of the ionic abundance Cl"^ , Cl^"^ and Cl^"^ . For 
argon, we have determinations of Ar'^^ and Ar''"'" but some 
contribution of Ar^ is also expected. We have adopted the 
ICF recommended bv llzotov et al.l l| 19941 ) for this element. 
Finally, we have used an ICF scheme based on photoion- 
ization models of [Rodriguez fc RubinI ((20051) to obtain the 
total Fe/H ratio. In Table[3]we include the total abundances 
determined for t^ = and t^ = 0.045 ± 0.007. The varia- 
tions due to the dependence of the adopted ICFs on the t^ 
considered are also included in the total abundances given 
in Table El 

In Table [3] we also include the ionic and total abun- 
dances obtained bv lCopetti et al . ( 2007) for their ABC spec- 
trum of NGC 2579. Those authors do not consider tem- 
perature fluctuations. We can see that all abundance s of 
twice ionised species determined bv lCopetti et al.l (120071) are 
somewhat larger than our determinations. This fact is per- 
haps indicating that both sets of observations correspond to 
nebular areas with somewhat different integrated ionization 



conditions. However, the total abundances obtained in both 
studies are in general consistent except in the case of Ne and 
CI, where the difference is 0.33 and 0.34 dex, respectively. 
The discrepancy in Ne/H ratios is due to the large difference 
between the intensities of the [Ne in] lines reported in both 
datasets. In the case of Cl/H the dif ference could be due 
to the different methodology used by ICopetti et al.l (|2007l ) 
and us for determining the abundance of this element. Those 
authors use Cl^'''/H"'" and an ICF and we determine Cl/H 
directly from the sum of the different ionic abundances. 

It is interesting to compare the abundance ratios de- 
termined for NGC 2579 with those of reference objects. In 
Table [J we show the C/0, N/O, Ne/O, S/O, Cl/O and 
Ar/0 measured for NGC 2579 in the case of assuming t^ 
— and t ^ = 0.045, the abund ance ratios for the solar pho- 
tosphere (JAsplund et al.ll2009l ) - the protosolar values are 
identical - and the values obt ained for the Orion nebula 
(|Garcia-Roias fc Estebanll2007l . assuming t^ > 0). In gen- 
eral, we can see that the ratios are rather consistent in the 3 
objects within the errors. The value we obtain of log(Ne/0) 
~ —0.93 is somewhat low, even lower than the typical ratio 
of about —0.70 found for extra galactic H ll regio ns covering 
a very wide range of metalicity (|Dors et al.ll2013l ). This indi- 
cates that our determination of Ne^'''/H"'" and Ne/H is prob- 
ably about 0.20 lower than expected. In the c ase of the N/O 
ratio, it is consistent with that obtained by ICopetti et al.l 
((2003) within the errors, but about 0.17 dex lower than that 
of the Orion nebula. 



5 CHEMICAL EVOLUTION MODEL AND 
DISCUSSION 

The derivation of precise abundances for NGC 2579 is an op- 
portunity to explore the behaviour of the radial abundance 
gradients of ionised gas at galactocentric distances larger 
than the solar ones. In Figure Owe include the C/H, O/H 
and C/O r atios of NGC 2579 as tho s e for t he sample of Hll 
regions of ICarcia-Roias fc EstebanI (|2007|), which are also 
determined using the intensity of RLs and the same method- 
ology. The objects of that sample are located at galactocen- 
tric distances between 6.3 and 10.4 kpc (assuming the Sun at 
8 kpc) . We also include the radial g radients determined from 
the linear fit to the data sample of lCarcia-Roias fc EstebanI 
(|2007i ) but extrapolated up to 13 kpc. As it can be seen 
in Figure (3] the point of NGC 2579 stands above the ex- 
trapolated gradients, specially in the case of the C and C/O 
gradients. Assuming that the distance to this object is well 
established (that seems to be t he case, see discussion of sec- 
tion 5.2 of ICopetti et alj|2007l ). our data indicate that the 
chemical behaviour of NGC 2579 does not fit the trend of 
the inner Galactic Hll regions. Although we deal with data 
of a single object - and obviously lacking statistical signifi- 
cance - this result indicates that the chemical composition 
of NGC 2579 is consistent with flatten gradients at its galac- 
tocentric distance. Another interesting feature of Figure O is 
that the C/O gradient seems to be better represented by a 
step function instead of a linear fit. Regions with Rg < 7.3 
kpc show log(C/0) ratios of the ionised gas of about 0.0 but 
regions farther out from the Galactic Centre can be repre- 
sented by a constant value around —0.26 with a rather small 
statistical dispersion. 
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Table 4. Comparison of abundance ratios 





NGC 2579 




Orion ncb.*^ 




(2 = 


i2 >0 


Sun^' 


i2 >0 


c/o 


— 


-0.26±0.06 


-0.26±0.07 


--0.28±0.04 


N/0 


-1.00±0.04 


-1.01±0.06 


-0.86±0.07 


~0.84±0.10 


Ne/O 


-0.96±0.02 


-0.93±0.08 


-0.76±0.11 


-0.77±0.08 


S/O 


-1.67±0.02 


-1.62±0.07 


-1.57±0.06 


-1.48±0.05 


Cl/O 


-3.47±0.02 


-3.47±0.06 


-3.19±0.33 


-3.42±0.05 


Ar/O 


-2.14±0.02 


-2.17±0.05 


-2.29±0.14 


-2.08±0.07 



■^ Asplund et al. (2009) . 

*> Garcia-Rojas & Esteban (2007). 




of slope of the gradients at large distances. Strong evidence 
for flat radial gradients from the spectra of Hll regions 
in the onter discs of nearby spiral galaxies has been pre- 
sente d in several recent pap ers by Brosolin ot al. (2 009) for 
M83 iGoddard et al.1 (|201ll ) for NGC 4 625. iBresolirTet all 
(|2012l ) for NGC 1512 and NGC 3621 and lWerk et al.1 (|201ll ) 
for a sample of 13 - mostly interacting - galaxies. More- 
over, a flattening or an upturn of the metallicity gradient 
has also been found from stellar photometry of red giant 
branch stars in th e outer disc o f the s p iral g alaxies NGC 300 
and NGC 7793 (|Vlaiic et al.l |2009| . I2OI1I ). The value of 
12-|-log(0/H) - obtained from the intensity of CELs - in 
the outer discs of the above mentioned galaxies ranges be- 
tween 8.2 and 8.4 (Bresolin et al. 2012 ), cons istent with the 
value of 8.31 ± 002 we derive for NGC 2579. iBresohn et aD 
(l2012f ) indicate that the flattening of the radial abundance in 
external spiral galaxies occurs approximately at the isopho- 
tal radius, -R25. In the case of the Milky Way, R25 is dif- 
flcult to estima te. The available determinatio ns give val- 
ues of 11.5 kpc ( de Vaucouleurs fc Pencelll978l ) or 13.4 kpc 
(IGoodwin et alJ 1998h. close t o the galactocentric distance 



determined by ICopetti etaP (J2007l l for NGC 2579. Metal- 
licity studies based on Cepheids and open clusters indicate 
that the cha nge of the slope o f the Fe/H ratio occurs at 
about 9 kpc (JLepine et al.ll201ll ). 

We have made a first exploration of the possible reasons 
of the flattening of the abundance gradients in the outer 
galactic discs buildin g a tailored chemical evolution model. 
iBresolin et al.l (J2012l ) discuss several mechanisms that can 
produce such phenomenon: flattening of the star formation 
efficiency, radial metal mixing or enriched infall. In this 
paper, we focus our attention on the first possibility be- 
cause it is the simplest one to implement in our available 
models. The model we have used is almost identical to the 
int ermediate wind yie l ds ori e (IWY), explained in details 
by ICarigi fc PeimbertI |2oi3). That IWY model was built 
to match three observational constraints along the Galactc 
disc: the radial distributions of the surface density of the to- 
tal baryonic mass, Mtot{r), and gas mass, Afgas(»"), and the 
radial gradient of the 0/H ratio defined by Hll regions of 
the Milky Way located between 6.3 and 10.4 kpc and derived 
using RLs. 

The main assumptions and ingredients of the IWY 
model are described in the following. The Milky Way 
disc was formed in an inside-out scenario from primor- 

^ dial infall with time scales r = r(kpc) —2 Gyr. T he ini- 

Yong et al.l 12012! ) . provide further evidence of such change tial mass function is that by iKroupa et al.l (|l993l ) in the 



Galactocentric distance (kpc) 

Figure 3. Ionised gas phase O, C and C/O radial abundance gra- 
dients of the Galactic disc from H II region abundances determined 
from recombination lines. The point at 12.4 kpc corresponds 
to NG C 2579, the rest are taken from lCarcia-Roias &: EstebanI 
(|2007l V The lines indicate the radial gradients obtained from the 
linear fit of the data between 6.3 to 10.4 kpc and extrapolated to 
13 kpc. The Sun is located at 8 kpc. 



Some authors have claimed - based on both H 11 regions 
and planetary nebulae data - that gas phase radial abun- 
dance gradien ts may fiatteri out at the outer parts of the 
Galactic disc (iFich fc Silkevlll99ll:IVilchez fc Estebaiill 19961: 



Maciel fc Quirezal Il999l : ICosta et al.1 |2004| : iMaciel et all 
20061) . although other authors do not support s u ch flattening 
JDeharveng et al.1 I2OO0I: iRudolph et al.1 I2OO6I : iHenrv et al.] 



l20ld : iBalser et al.l 120111 ). Metallicitv 



other kinds of ob j ects, s uch as Ceph e ids (iLtick et al 



determinatio ns of 
20031: 



Andrievskv et al.1 |2004| : lYong et al.1 I2OO6I: iPedicelh et al.1 



2OO9I) and Gala ctic open clusters ( Twarog et al.l 119971 : 
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0.08 — 80 Mq range. We consider an array of metal de- 
pendent yields: (a) for low and inter mediate mass stars 
(0.8 ^ m/ Mp, < 8), we use the yi elds bv lMarigo et all l|l996l. 
Il998l ) and IPortinari et al.l l|l998l ): (b) for bin ary stars (3 ^ 
mbin/MQ ^ 16), we adopted t he yields by Nomoto et al.l 
lH^) in SNIa formulation by iGreggio fc Renzinil (|l983l ). 
We used Aun = 0.08, as the fraction of binary stars 
that are progenitors of SNIa; (c) for massive stars (8 ^ 
m/MQ ^ 8 0) we conside red t he yields including stellar 
rotation by iHirschil (|2007l ') and iMevnet fc Maedeil (|2002l ') 
for Z ^ 0.004; and the interme diate wind yiel ds, obtained 
as an average of the yields by J Maederl (|l992l . high mass- 
loss rate) and iHirschi et al.l l|2005l . low mass- loss rate)) for 
Z = 0.02. Since Fe yields are not computed by the Geneva 
gro up, we adop ted Wooslcy & Weaver ( 1995*) values, follow- 
ing [Caxigi_&_jjernande^ (|2008l ) prescription. 

Our model differs from IWY only in the assumption of 
a variable star formation process efficiency, i^, for the outer 
galactocentric radius. With this parameter, the star forma- 
tion rate, SFR, is a spatial and temporal function of the 
form 

SFR(r,t) = u{r) x M^;,t{r,t) x (Mgas + M,taO°-''(r, t); (2) 

where Mstar is the surface mass density of stars. The for- 
mation process efficiency f was kept constant in the IWY 
model, while in the present paper we assume that u de- 
pends on the galactocentric radius. This is the only differ- 
ence between the IWY model and the present one. We adopt 
u = 0.016 = i^in for r < 10 kpc and u/uin = 0.8, 1.3 and 2.2 
for r =10, 12 and 14 kpc, respectively. The values of u as 
a function of the galactocentric radius we have chosen are 
those that reproduce the behaviour of O/H radial gradient 
at large galactocentric radii. 

In Figure [l] the model is compared to observational 
constraints, as follows: 

a) Mgas{r) values are obtaine d by adding the atomic an d 
molecular data from figure 7 bv lKennicutt fc Evand l|2012h . 
These original data include the hydrogen and helium com- 
ponents. 

b) SFR(r) values are t aken from the same figure 7 by 
iKennicutt fc Evans! (|2012l '). 

c) The disc scale lengt h, Rd, of t he Ga lactic stellar disc 
{Mstar{r) oc e-''/^'^) bv lYin et all l|2009l ). In their table 7 
they show Rd values between 2.3 and 5 kpc for K, R, V and 
B bands. 

d) O/H and C/H values from Hll regions and com- 
puted from RLs IjGarcia-Roias fc Estebarj|2007l . this work). 
The total gaseous abundances shown in our Table [3] have 
been modified under the assumption that 35% of the O 
atoms and 25% of the C at oms are trapped in dust grains 
l|Peimbert fc Peimbertll2010l ). increasing the O/H and C/H 
values by 0.12 and 0.10 dex, respectively. 

e ) N/H values computed from CELs and assuming t^ > 



i Garci'a-Roias fc Estebanli2007 . this work). We have as- 
sumed no dust depletion for N. 

f) Fe/H grad i ents fr om open cl u sters and Cepheids by 
iPediceUi et al] l|2009l ). iLuck et all (|201ll ) and lYong et al.1 
1 20121) . 



solute values) and with the SFR, Mgas and Mstar radial 
distributions, for which the model was not built at all. 

The results of the chemical evolution model we have 
built in this work indicate that an increase of v for r ^ 10 
kpc can explain the chemical flattening in the outer Galac- 
tic disc. It is important to remark that v does not represent 
the star formation efficiency, SFE, defined as SFE(r) = 
SFR{r)/Mgas{r), the SFR per unit of surface mass density 
of gas. In order to clarify this point, we show in Figure [5] 
the radial distribution of v and SFE. In this figure, we note 
that a constant value of f at inner radii implies a decrease of 
SFE and an increase of f at outer radii implies a flat SFE. 
A sim ilar behaviour of the SFE{r) was found bv lBigiel et al.l 
(2010!), based on the combination of atomic hydrogen and 
far-ultraviolet emission in 17 spiral and 5 dwarf galaxies. 
They conclude that the SFE in the outer disc (r > -R25) is 
flatter compared to the SFE in the inner disc (r < -R25), 
due to the increase of the SFR with H I column density. The 
SFE is the inverse of the gas depletion time, the time re- 
quiered for present-day star formation to consume the avail- 
able gas. The values of SFE we obtain at the outer disc 
indicate a depletion time of about 7.5 Gyr and that this 
value becomes alm ost constant in that part of the Galaxy. 
iBigiel et al.l (|2010|) find depletion times between 10 to 100 
Gyr for their sample of spira l and dwarf galaxies. 

As lBresolin et al.l (|2012l ) point out, a flat SFE can ap- 
proximately translate into a flattening of the abundance gra- 
dient, using their equation 1: 



O 

H 



{yoxtx Ssfr) 
fi X Shi 



(xSFE; 



(3) 



It is remarkable the agreement with the observational 
constraints for r ^ 10 kpc, specifically with the current C/H, 
N/H, C/O and N/O abundance gradients (slopes and ab- 



where yo is the net oxygen by m ass (0.006, integrated yiel d 
obtained from the IWY model, ICarigi fc PeimbertllJOllI '). 
fi = 11.81 corresponds to the conversion factor from abun- 
dance by number to abundance by mass; t is the timescale for 
the star formation and "Esfr and T.hi are the SFR and H I 
surface density, respectively. Assuming a constant SFR and 
gas density and that they have been equal to the present-day 
values in the last Gyrs, we can estimate t from the equation 
given above. From Figure|31 at the distance of NGC 2579, we 
obtain T.sfr ^ 1 Mq pc"^ Gyr"^ and Eh/ -- T.gas ~ 8 Mq 
pc~^. With these numbers we obtain that the time required 
to enrich the ISM up to the oxygen abundance of NGC 2579 
determined from CELs and adding 0.12 dex due to dust de- 
pletion, 12 + log (0/H)cELs = 8.43, is ~ 4.2 Gyr (3.2 Gyr 
without dust correction). In the case we use the abundance 
determination based on RLs including dust correction, 12 -|- 
log (0/H)rls = 8.67, the time is - 7.5 Gyr (5.7 Gyr with- 
out dust correction) . Simulations of galaxy formation in a A 
cold dark matter (ACDM) universe predict that the outer 
discs of galaxies could have mean formation times around 
4-6 Gyr, while the inner parts should have mu ch longer val- 
ues, > 10 Gyr (|Scannapieco et al.ll20oi , l2009l ). Our estima- 
tions for the outer Galactic disc a re roughly consistent with 
these predictions. iBresolin et al.l (2012) determine enrich- 
ment timescales for the external discs of two spiral galaxies, 
for the apparently isolated spiral NGC 3621 they obtain a 
time of '^ 10 Gyr, and a comparatively much shorter value - 
around 2-3 Gyr - for the interacting spiral galaxy NGC 1512. 
In an upcoming paper, we explore other factors that 
may produce flat gradients in the outer Galactic disc, as ho- 
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Figure 4. The present-day radial distribution of surface mass densities of gas (Mgas) and stars (Mstar) in Mq pc~ units and star 
formation rate (SFR) in Mq pc~^ Gyr~^ units (left panels); ISM abundance ratios and stellar Fe/H gradient (right panels). Continuous 
black lines; Results of the chemical evolution model that assumes an increase of u for r ^10 kpc. The three dotted lines parallel to the 
predicted Mstar correspond to scale lengths (Rd) of 2.8, 2.3 y 3.4 kp c for 4-8 kpc, 8-10 kpc and 10- 14 kpc, respectively. Observational 
data: filled circles represent Hll regions (corrected for dust depletion. iGarcia-Roias fc E steban 2007, this work); red circles: NGC 2579 
(corrected f or du st depletion, this paper); dis continuous line s : Fe/ H gradient by lYong et al. (2012, short-long-dashed magenta line), 
iLuck et al.l ll201ll . short-dashed blue lin e) an d iPedicelli et al.l 1 120091 . long-dashed green line). Filled green squares: values of Mgas and 
SFR compiled by|Kcnnicutt & Evans (2 0121 ). Vertical bars: average observational errors. Diagrams showing O/H, C/H, N/H and Fe/H 
express the abundances in units of 12-|-log(Xi/H) and those showing C/O and N/O are in log(Xi/0) units. 



mogeneous infall, enr iched accret ion, and gaseous and stellar 
radial migration f see iBresolin et al..,201Z . for references). 



6 CONCLUSIONS 

We present deep echelle spectrophotometry in the 3550- 
10400 A range of the Galactic H ll region NGC 2579. This 
object is located at a rather large g alactocentric distance, 
12.4 ± 0.7 kpc (jCopetti et al.l |2007| ) and has been largely 
neglected due to ide ntification problems finally resolved by 
ICopetti et al.l (|2007l ). It also shows a rather high surface 
brightness and ionization degree. All these properties make 
NGC 2579 an excellent object to explore the behaviour of 



the Galactic radial abundance gradients of the ionised gas 
in the outer disc of the Milky Way. 

We have derived consistent and precise values of the 
physical conditions of the nebula making use of several emis- 
sion line-intensity ratios as well as abundances for several 
ionic species from the intensity of coUisionally excited lines 
(CELs). Our deep spectra permits - for the first time for 
a Galactic Hll region with so low metallicity - to obtain 
a very good determination of the ionic C^"'', O^ and O^^ 
abundances - as well as the total O abundance - from the 
intensity of faint pure recombination lines (RLs) . The com- 
parison between the 0^/H^ and 0^^/H"'" ratios determined 
from RLs and CELs gives an abundance discrepancy factor 
(ADF) of 0.15 ± 0.08 and 0.27 ± 0.03 for 0+ and 0^+, re- 
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Figure 5. The present-day radial distribution of: v, where 
u{r) = SFR{r) X Mga^ir) X (Mgas + M star)-° ■* (r) , in 0.01 
Gyr~^(M0 pc~'^)~^'^ units (upper panel), the star formation 
efficiency, SFE{r) = SFR{r)/Mgas{r), in Gyr"! units (middle 
panel) and depletion time TDcp('') = 1/ SFE(r) in Gyr (lower 
panel). 



spectively. Values in complete agreement to the ADFs found 
for other Galactic and extragalactic Hll regions. The value 
of t^ that produces the agreement between the O^"'' abun- 
dance determined from CELs and RLs is t^ — 0.045 ± 0.007, 
which is consistent to the values we obtain from the com- 
parison of the O^ abundances - t^ — 0.048 ± 0.029 - and 
applying a maximum likelihood method for minimising the 
dispersion of the He^/H"*" ratio from individual lines. 

Our abundance results for NGC 2579 based on RLs 
permit to extend the previous determinations of the C, O 
and C/O gas p hase radial gradients of the inner Galactic 
disc obtained bv lEsteban et al.l (|2005h to larger galactocen- 
tric distances. We find that the chemical composition of 
NGC 2579 is consistent with flatten gradients at its galac- 
tocentric distance. This result is in agreement with previous 
claims of flatten outer abundance gradients based on results 
for ionised nebulae, Cepheids, and open clusters in the Milky 
Way and with recent results for the outer discs of nearby 
spiral galaxies. 

We have built a chemical evolution model that repro- 
duces the observed radial distributions of O, C and N abun- 
dances, the C/O, N/O ratios and other observational con- 
straints, finding that the flattening of abundance gradients 
in the Milky Way outside the isophotal radius can be ex- 
plained due to a shallowed star formation efhciency for 
r ^ 10 kpc. 
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